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Abstract. The contribution of the different terms of the running-coupling Balistky- 
Kovchegov (rcBK) equation to the description of inclusive HERA data is discussed. 
Within this framework an alternative definition of the saturation scale is presented. 
The definition is based on the ratio of the linear to the non-linear term of the 
equation. A similar ratio is used to study the contribution of the non-linear term 
to the evolution of inclusive HERA data with rapidity. It is found that, although the 
data are well described, the behaviour of the different terms of the rcBK equation for 
HERA kinematics is not what it is naively expected from saturation arguments. 


1. Introduction 

One of the most interesting problems nowadays in QCD is to understand the high-energy 
behaviour of the gluon structure of the proton. Inclusive data from deeply inelastic 
scattering at HERA show that the proton structure function F 2 {x,Q^) grows rapidly 
with decreasing values of Bjorken-x at hxed values of the virtuality, Q^, of the photon 
[1]. At small values of x, which correspond to the high-energy limit of QCD, the main 
contribution to F 2 {x, Q^) comes from the gluon density. The growth of the gluon density 
for decreasing x at a fixed scale was predicted in [a El SI El E], where it was understood to 
be due to gluon-branching processes. The unitarity of the cross section implies that the 
gluon density should stop growing at some point. This is known as saturation. It was 
predicted to occur in QCD in the seminal work [7], where a nonlinear term was added 
to the evolution equation of the gluon density. The presence of saturation generates a 
hard scale, the so called saturation scale, Qs{x), which allows to perform perturbative 
calculations, in spite of the high density of the system. 

A framework to study saturation phenomena is the Color Glass Condansate (CGC); 
an effective theory, which describes the high energy limit of QCD Pi. The CGC 
generates a set of equations known as the JIMWLK equations [ini El Ea na El E]. 
These equations are equivalent to a hierarchy of equations found by Balitsky [ISIEI 
[Tallin]. In the limit of a large number of colours, the hierarchy reduces to one equation, 
which was independently derived by Kovchegov [201 E] within the colour dipole model 
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[221 Eai Eg. This equation, called the Balitsky-Kovchegov (BK) equation, is simpler 
than the JIMWLK equations, but it already embodies most of the relevant behaviour as 
shown by a comparison of numerical solutions in both approaches [25l ESj- Corrections 
to include the running of the coupling, were derived in EB] and the resulting 
equation, called rcBK, was solved numerically in [29]. The solutions of the rcBK 
equation have been used to describe a variety of HERA data: inclusive [30l |3T1 l32] . 
diffractive [32] and exclusive production of vector mesons [STl 133] . Furthermore, the 
rcBK equation includes naturally geometric scaling, one of the most striking discoveries 
with HERA data [3S]. Due to these successes the rcBK equation is one of the main 
tools used today to investigate the phenomenology of saturation in QCD. 

Here, we use the rcBK equation to study the rapidity dependence of saturation 
in inclusive HERA data. In section El we introduce the rcBK equation, review the 
numerical method to solve it and describe the behaviour of its solution, the dipole 
scattering amplitude N{r,Y), as a function of rapidity Y and transverse size r of the 
dipole. We also discuss an alternative definition of the saturation scale and compare its 
evolution in rapidity to that of the usual definition. In section [3] we present F 2 {x,Q‘^) 
in the colour dipole formalism and compare it to HERA data. We also analyse the 
evolution in rapidity of the contribution of the non-linear term to the description of 
F2 {x,Q‘^). In section m we comment on the geometric scaling properties of the dipole 
scattering amplitude and of the inclusive cross section for the scattering of a virtual 
photon off a proton. In section El we discuss the implications of our results. Finally, in 
section El we summarise our findings. 

2. The dipole scattering amplitude in the rcBK framework 

The BK equation describes the evolution in rapidity Y of the scattering amplitude 
N{r, Y) for the scattering of a colour dipole of transverse size r with a target. During 
the evolution the dipole may split into two dipoles or two dipoles may recombine into 
one. Here we follow the work of [211 EDI EH and study the equation without impact 
parameter dependance. 

2.1. The rcBK equation 

The running coupling Balitsky-Kovchegov equation [2Z112SI is 

dN^^Y) ^ I d^^p^^f^^^^^^^(^p^^^^^Y) + N{r2,Y)-N{r,Y)-N{ruY)N{r2,Y)j ( 1 ) 

where r 2 = r — ri; r = |r| and similarly for ri and r 2 . The kernel incorporating the 
running of the coupling is given by 

as{r^)Nc ( r^ 1 f \ 1 /a.(ri) \\ 

27r rf )^rj[aM) ))^ 


K{r,ri,r2) 


( 2 ) 
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with 


«s(r 2 ) = 


Att 


( 11 - 


iNf) In 


4C2 
r*2 A ^ 

^^QCD 


(3) 


where Nf is the number of active flavours and C is a parameter to be fixed by comparing 
to data. 

The rcBK equation can be separated into the linear, 

/^(r, Y)=J dnK{r, n, r^) F) + iV(r 2 , F) - N{r, Y)j , (4) 

and the non-linear part, 

lNL{r,Y) = J dnK{r,f,,f2)N{ruY)N{r,,Y). (5) 

When the contribution oi Inl is small with respect to that of II, the evolution in rapidity 
of the dipole scattering amplitude is driven by the linear term. When both terms are of 
the same size, and equilibrium is reached, there is no more evolution and the scattering 
amplitude is said to be saturated. 

The final ingredient to the rcBK equation is the initial condition. For the initial 
form of the dipole scattering amplitude we use the McLerran-Venugopalan model 133: 

iV(r,y = 0) = l-exp(-F^ln(^ + e) 

with the values of the parameters QIq, C and 7 taken from fit (e) in Table 1 of [3T] . 
Note that in the fit the initial rapidity, F = 0, is at xq = 0.01, where the relation 
between F and a: is F = ln(a:o/a:). Within this prescription rapidity F = 7 corresponds 
to X ~ 10“®; that is, to the smallest x measured at HERA for perturbative scales. The 
fit yielded the following values for the parameters: = 0.165 GeV^, 7 = 1.135 and 

C = 2.52. The fit was performed under the assumption that as{r^) freezes for values of 
r larger than tq defined by as(rQ) = 0.7. We set Aqcd to 241 MeV. 



2.2. The dipole scattering amplitude 


The rcBK equation can be solved numerically using Runge-Kutta (RK) methods in 
parallel over a grid in r. This approach has been followed in [2911301 [3T] using a second 
order algorithm. We have solved the rcBK equation using the following RK methods: 
forward Euler (order one), Heun (order two) and Classical (order four). In all cases we 
have a grid of 800 points in r, spaced uniformly in ln(r) from r = 10“® to 100 GeV“^. 
The equation was solved using steps of 0.05 rapidity units. The formulae for the RK 
methods are explicitly given in [Appendix A 


The upper panel of Figure [T] shows the solution of the rcBK equation as a function 
of r at different rapidities for the RK methods of order one and four. (Heun’s method 
yields results indistinguishable from those of the Classical method, thus it has not been 
included in the figure.) The forward Euler method lags slightly behind the Classical 
method. It has been checked that our results agree with those of [31], where a second 
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order RK method has been used. All following results are based on the Classical method. 
The lower panel of Figure [T] shows the dipole scattering amplitude as a function of 
rapidity at several fixed dipole sizes. The dipole scattering amplitude at a fixed r shows 
a fairly constant behaviour from Y = 0 up to a rapidity which depends on the value 
of r, and from there onwards it raises steeply with Y until it saturates. For values of 
r smaller than around 0.1 GeV“^, the growth of A^(r = const., Y) with rapidity starts 
beyond the kinematic reach of HERA. Even for r = 5 GeV“^, N{r = 5 GeV“^,R) is 
fairly constant down to x ~ 5 ■ 10“^. 

2.3. The saturation scale 

In order to study the relative importance of the linear to the non-linear contributions 
to the evolution of the dipole scattering amplitude with rapidity we define the ratio 



(7) 


The upper panel of Figure [2] shows Rnl 2 L as a function of r for three different values 
of rapidity. In all cases, the ratio rises with increasing r. At higher rapidities there is a 
similar behaviour of the ratio and the dipole scattering amplitude: The shape of Rnl 2 L 
is the same at different rapidities, but it is shifted towards smaller values of r, just like 
in the upper panel of Figure [H The lower panel of Figure [2] shows Rnl 2 L as a function 
of rapidity at fixed values of the dipole size r. Again, the behaviour is very similar to 
that of the dipole scattering amplitude. The ratio is flat or even slightly decreasing up 
to a value of rapidity, which depends on r; from there onwards the ratio rises steeply 
until it reaches unity. Smaller r have a small ratio at the initial rapidity, while large r 
has a large ratio. 

One commonly used definition of the saturation scale is A^(r = 1/Qs, Y) = k, where 
K is an arbitrary value close to one. Given the behaviour of the ratio Rnl 2 L 5 there is 
another possibility to define Qs{x), which is perhaps more natural in the context of the 
rcBK equation: the saturation scale is the inverse of the dipole size, where the ratio 
of the non-linear to the linear contribution to the rapidity evolution of the scattering 
amplitude is k: 


RnL2l{'>^ — ^/Qs^Y) — 


( 8 ) 


Figure [3] compares the evolution with rapidity of the saturation scale for the two 
definitions. The qualitative behaviour is similar in both cases: the saturation scale raises 
very slowly for a few units in rapidity, then there is a large region where Qg changes 
faster until it reaches a power-law increase with rapidity, which sets in for Y k, 20-30. 
Quantitatively, at small rapidities, the saturation scale using equation ([7]) with k = 0.5 
is very close to the standard definition of the saturation scale with k = 0.25, while at 
large rapidities is closer to the standard definition with k = 0.5. 

It is clear then that both definitions are qualitatively equivalent. The commonly 
used definition is more general, while the one we introduce here is specific to equation 



N(r=const.,Y) N(r,Y=const.) 
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r(1/GeV) 



Figure 1. The solution of the rcBK equation is shown. Upper panel: dipole scattering 
amplitude N{r, Y) as a function of r at rapidities U = 0, 7, 50 and 100 (from right 
to left) computed with the forward Euler and the Classical methods. Lower panel: 
dipole scattering amplitude N{r,Y), found with the Classical method, as a function 
of rapidity at several values of r. 











RNL2L(''=const.,Y) RNL2L(''>Y=const.) 
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Figure 2. The relative contribution of the non-linear to the linear term of the rcBK 
equation to the evolution of the dipole scattering amplitude with rapidity is shown as 
a function of r at constant rapidity (upper panel) and as a function of Y at constant 
r (lower panel). 
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Figure 3. Evolution with rapidity Y of the saturation scale (5s(h") for the two 
definitions discussed in the text. 

©■ As this latter definition is tailored to the rcBK equation, it is helpful to understand 
qualitatively some of the properties of equation ([T]). 

3. Description of inclusive HERA data in the rcBK framework 

3.1. The proton structure function F 2 {x, Q^) in the colour dipole formalism 

The proton structure function F 2 is related to the cross section a^P for the scattering 
of a virtual photon 7 * off a proton p by: 

F2{x, Q^) = Q^) + , (9) 

where the subindices T and L refer to the transverse, respectively longitudinal, 
polarisations of the virtual photon. 

In the colour dipole picture [22l |23l |2l] , the incoming photon splits into a quark- 
antiquark pair, which forms a colour dipole, and this dipole interacts with the target. 
For the case of no impact-parameter dependence the cross section is given by 

afy(x,Q^) = o-Qy drj^ dz\^iy'^\z,r,Q‘^)fN{r,x). (10) 
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To lowest order in aem the square of the wave functions describing the splitting of the 
virtual photon into the colour dipole are given by 

and 

- ^)^Al(£r)), (12) 

where we have follow the normalisation convention of [38]. In these equations, 2 ; is the 
momentum fraction of the photon carried by the quark, Kq i are Macdonald functions, 
6 / is the electric charge of the quark of flavour / in units of the electron charge and 

= z{l — z)Q‘^ + mj. (13) 


Only three quark flavours have been considered. For each flavour, the parameter my has 
been hxed at 140 MeV/c^ and the value of (Tq = 32.895 mb has been obtained by the £t 
mentioned above. Finally, following the convention introduced in [32], which allows to 
have a smoother approach to the photoproduction limit, we set 


X = x\ 1 + 


Amj 


(14) 


The upper panel of Figure H] shows 

|vl/(r,g2)|2^ + (15) 

for four different values of the virtuality of the photon, Q^. The wave function grows 
rapidly for decreasing dipole sizes reaching the same value for all shown values of at 
r 0.1 GeV“^. For large dipoles there are sizeable differences among the wave functions 
corresponding to the different virtualities. The lower panel of the same hgure shows the 
integrand of equation ffTOjl at T = 2. The factor originates in the Jacobian to go 
from df to dr and from dr to (Jln(r). As expected, the integrand goes to zero at small 
r approximately as r^. The position and height of the maximum of the curves depends 
on Q^: for the smaller value of the wave function |Tp peaks at larger dipole sizes 
than for larger values. For these larger values of the shape of |Tp is less peaked, 
more broad. Qualitatively, the same behaviour is found at larger rapidites. 


3.2. Rapidity evolution of the non-linear contribution to F 2 {x,Q^) 

The upper panel of Figure [5] shows the structure function F 2 {x,Q‘^) as measured by 
HERA [1] for four different values of the photon virtuality along with its description 
using the colour dipole formalism and the solution of the rcBK equation. The relation 
between rapidity and x is given by F = ln(a:o/x), where xq = 0.01 and x as been dehned 
in equation fll4p . Note that the data at = 200 GeV^ were not included in the fit, 
so in a sense the corresponding curve is predicted by the formalism. The description of 





r \W{r,Q )| N(r,Y=2) \W{r,Q )| (GeV ) 
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r(1/GeV) 

Figure 4. Upper panel: The sum of the transversal and longitudinal wave functions 
integrated over z, see equation (I15I) . is shown for four different values of the virtuality 
of the photon, Q^. Lower panel: The integrand of equation (fTUl) at Y = 2 is shown for 
different values of Q^. 
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data by the model is quite good over three orders of magnitude in and several units 
of rapidity. 

Looking at equations ([9]) and ffTOl) it is clear that the evolution of F 2 with rapidity 
is driven by the evolution of the dipole scattering amplitude. The left-hand side of 
equation ([1]) allows to separate the contribution to the evolution of F 2 coming from the 
linear and the non-linear terms. We dehne dFf^^/dY and dFf/dY as the result of 
exchanging N{r,Y) by lNL{r,Y), respectively Iiir^Y), in equation fITUD . Their ratio is 
shown in the lower panel of Figure [5l The ratio decreases slightly with rapidity for a 
few units, and then starts to grow with rapidity towards one. To the lowest values of 
corresponds the larger values of the ratio: above 0.5 for = 0.2 and 2 GeV^ and 
below 0.5 for larger values of Q^. 

4 . Geometric scaling 

One of the most striking phenomena discovered with HERA data is geometric scaling 
[36] . On the other hand, the rcBK equation includes naturally geometric scaling in 
its solutions at large rapidities. It has to be noted that these two geometric scalings 
are not the same, that is the scaling variables are different, but are intimately related. 
The upper panel of Figure [6] shows the dipole scattering amplitude obtained solving the 
rcBK equation as a function of r = rQs(Y) at different rapidities. For the saturation 
scale, the dehnition given in equation ([8]) has been used with n = 0.5. In this case, 
geometric scaling means that the dipole scattering amplitude at different rapidities is 
the same when plotted as a function of r. It can be seen that the dipole scattering 
amplitude at the initial condition is far away from this behaviour. Even at K = 10 
the dipole scattering amplitude does not yet display geometric scaling. Only for larger 
values of rapidity the geometric scaling regime is approached. The lower panel of the 
same hgure shows the behaviour of the 'y*p cross section as a function of the scaling 
variable Q'^/Q^fY). For = 0.2 GeV^, Q'^/Q^fY) is approximately 0.3 at K = 7 and 
grows by a factor of ten by each order of magnitude increase in Q^. The different curves 
collapse into each other only for values of Q‘^/Q‘l{Y) less than 0.01, that is for rapidities 
larger than 7. 

5. Discussion 

The dipole scattering amplitude found with the rcBK equation gives a nice description of 
HERA data, which has been taken as supporting evidence of saturation effects already 
for HERA kinematic region. This argument is very strong. Geometric scaling is a 
striking behaviour, which is naturally found in saturation models. On the other hand 
the behaviour of a solution of a differential equation has two components, the dynamics 
inherent to the equation and the influence of the initial condition. 

The rise of the dipole scattering amplitude with rapidity for a hxed dipole size 
and the appearance of geometric scaling are behaviours associated with the saturation 


3yF2HF2 F (X,Q ) 


Rapidity dependence of saturation in inclusive HERA data with the rcBK equation 11 




Figure 5. Upper panel: description of HERA data from [T] using the colour dipole 
formalism and the solution of the rcBK equation. Here, Y — ln(0.01/i) where x as been 
defined in equation m- Lower panel: relative contribution of the non-linear to the 
linear term of the rcBK equation to the evolution with rapidity of F 2 for four different 
values of the photon virtuality, Q^. Here we have used the notation dy = d/dY. 











(Q /Q (Y)) (1/GeV ) 
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q"/q"(Y) 


Figure 6. Upper panel: The dipole scattering amplitude is shown as a function of 
T = rQs(Y) at different rapidities. For the saturation scale, the definition given in 
equation ([5]) has been used with k = 0.5. Lower panel: 7 *p cross section as a function 
of the scaling variable Q‘^/Ql{Y) for different values of Q^. 
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regime. These behaviours appear only at rapidities beyond those reached at HERA 
when using, in the initial condition, the parameters which allow a good description 
of data. This behaviour is not found only with this set of parameters. For example, 
[2U] reports similar observations when using very different set of parameters for the 
initial condition. The saturation scale behaves similarly. The behaviour expected in the 
saturation regime is reached only at large rapidites. 

Turning to the description of data, it is true that the fit is good. But there are two 
points to consider. First, the fit used data whose rapidity was already very close to the 
initial condition. This could have given an extra weight to the importance of the initial 
conditions with respect to the dynamical component of the rcBK equation. Second, 
data with very low values of have been included in the fit. (The fit included all 
available data with 0.045 < < 50 GeV^.) It is not clear if a perturbative approach is 

valid at these low scales. It is also not clear how important are these data to determine 
the values found by the fit. It can be seen in the lower panel of Figure [5] that it is 
only for these data at low that the non-linear term is numerically dominant for the 
evolution with rapidity, while at larger values of the evolution is driven mainly by 
the linear process. Another observation is that the linear process gains in importance 
with rapidity over the range where data exist, which is the opposite behaviour of what 
is naively expected of the dynamics of the rcBK equation. 

6. Summary and conclusions 

We have studied the contribution of the different terms in the rcBK equation to the 
description of inclusive HERA data. We find that although the data is well described, 
the behaviour of the model in the kinematic range of the measurements is not what is 
expected from the dynamics of the rcBK equation. One could try new fits to HERA 
data without inclusion of the measurements at the lowest values, which may not be 
adequate for a perturbative model, and give some time for the evolution of the solution 
before the first data to be fitted in order to reduce the weight of the initial conditions 
and to increase the importance of the dynamic content of the rcBK equation. 

These results point also to the need of having data at larger rapidities and at 
perturbative scales. In this context, the measurement of exclusive photoproduction of 
J/'^ at LHC both off protons [39l 001 01] as off nuclei 021031011 at higher energies 
than those reached at HERA is very good news. Ultimately the realisation of dedicated 
electron-hadron colliders 0306107] would provide high quality data to continue these 
studies and conclude if a saturation regime has been reached. 
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Appendix A. Runge—Kutta methods to solve the rcBK equation 

Using the following definitions 

Io = J dm{r,n,r2)iyN{r^,Y) + N{r2,Y) - N{t,Y) - N{r^,Y)N{r2,Y)^, (A.l) 

h= J df,K{f,f,,f2)(^l-N{r,,Y)-N{r2,Y)y (A.2) 

h = J dfiK{r,fi,f2), (A.3) 

the change of N{r,Y) after a step in rapidity AY can be found using Runge-Kutta 
methods over a grid in r. 

For the forward Euler method, the solution reads: 



A(r, Y + AY) = N{r, Y) + AYR. 



(A.4) 

For Heun’s 

method, the solution reads: 





A(r, Y + AY) = A(r, Y) + AYR + 

r r 

2 

{m\ 2 r 

2 

(A.5) 

For the Classical method, the solution reads: 





AY f 

N{r,Y + AY) = N{r,Y)+ ^ (Jo 

+ 2K2 + 27^3 

+ -^ 4 ^ , 

(A, 6 ) 

where 






AU (AU)2 „ 

K 2 ^R + -^Rh - 



(A.7) 


AY (AU)2 „ 

Ks — R + 2 -^ 27 i ^ -^2/2, 



(A. 8 ) 

and 






A4 = Jo + AUA 3/1 + {AY)^KlR. 



(A.9) 


In our implementation, the values of the dipole scattering amplitude for sizes r, which 
are not in the grid, are found by a linear interpolation in ln(r) between adjacent values of 
r in the grid. For dipole sizes larger than the larger size in the grid the dipole scattering 
amplitude has been taken as one, and for dipole sizes smaller than the smaller size in 
the grid the dipole scattering amplitude has been linearly scaled in r from the value 
corresponding to the smallest size in the grid. 
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